The Domino Effect in Decentralized Wireless 

Networks 



Brage Ellingsaster* and Torleiv Masengt 
*UNIK- University Graduate Center at Kjeller, University of Oslo, Norway 
t Norwegian Defence Research Establishment 

brage@unik . no 



Abstract — Convergence of resource allocation algorithms is 
well covered in the literature as convergence to a steady state is 
important due to stability and performance. However, research 
is lacking when it comes to the propagation of change that occur 
in a network due to new nodes arriving or old nodes leaving or 
updating their allocation. As change can propagate through the 
network in a manner similar to how domino pieces falls, we call 
this propagation of change the domino effect. In this paper we 
investigate how change at one node can affect other nodes for 
a simple power control algorithm. We provide analytical results 
from a deterministic network as well as a Poisson distributed 
network through percolation theory and provide simulation 
results that highlight some aspects of the domino effect. The 
difficulty of mitigating this domino effect lies in the fact that to 
avoid it, one needs to have a margin of tolerance for changes 
in the network. However, a high margin leads to poor system 
performance in a steady-state and therefore one has to consider 
a trade-off between performance and propagation of change. 

I. Introduction 

Convergence to a stable state is desirable when considering 
resource allocation in decentralized wireless networks since 
rapid transmitter configuration is a costly operation in terms 
of overhead and energy consumption. Stable states (or equi- 
libria) also allow us to analyze the performance of specific 
resource allocation algorithms efficiently in terms of the stable 
states reached by the algorithms. Thus much research has 
been devoted to the analysis of convergence of resource 
allocation algorithms 1 1 1|2||3]|4||5|. The desirable feature of 
such algorithms is fast convergence rate, convergence to an 
efficient stable state and convergence independent of the initial 
conditions. 

An aspect overlooked in the literature on physical layer 
resource allocation algorithms is the consequence of a change 
in the network and how it affects the allocation of other nodes 
in the network. Imagine a new node enters the network or an 
old node leaves or updates its allocation (e.g. changes power 
to satisfy a new SNR requirement), how does this affect the 
nodes in the network? Specifically we are interested in the 
propagation of change in networks due to such updates. When 
an update occurs, change propagates through the network in a 
way similar to a "domino" effect, changing the allocation of 
nodes in the network. 

It is desirable that resource allocation algorithms produce a 
minimal domino effect, meaning that as few nodes as possible 
must adapt to a change in the network. Considering large scale 
networks on the order of a hundred links, if an update at one 



link results in all other links changing their allocation as a 
consequence, stable states in the network will not occur if 
convergence time is longer than the frequency at which updates 
occur. One way to mitigate such a scenario would be for all 
nodes to allocate resources with a buffer or margin, similar to 
a fading margin. However, adding such a margin will lead to 
a loss in performance in stable states, and it is also difficult 
to say how large this margin should be, as it depends on a 
number of variables. Thus, we postulate that there will have 
to be a trade-off between the performance of steady states and 
how large a domino effect one can tolerate. To the best of 
the authors' knowledge, no work has been published in the 
literature which investigates these issues. 

We start the investigation by presenting a deterministic net- 
work where we can analytically characterize the domino effect. 
We then investigate a network where nodes are distributed in 
the plane according to a Poisson Point Process. One of our 
main results are obtained through percolation theory, which 
has been applied to wireless communication to study the 
connectivity of wireless multi-hop networks QQ. We show 
that there exists a critical density of nodes in the network so 
that for densities lower than this critical value, the propagation 
of change is finite almost surely (a.s.), and for densities larger 
than this critical value, the propagation of change is infinite 
a.s. Through simulation results we also note that the rounds 
in which nodes are affected follow a power law distribution 
for low margins. We can state that for these low margins the 
system is in a critical state so that the dynamics of the network 
is similar to systems exhibiting self-organized criticality (SOC) 
in physics [8]. 

A. Cellular vs. Decentralized 

It is important to understand why the domino effect is 
an important aspect when it comes to designing resource 
allocation algorithms for decentralized networks. We do this 
with a comparison to cellular networks, which do not have the 
problem of a domino effect. In a cellular network the resources 
are fixed and controlled strictly by the base station controller. 
E.g. in a CDMA network, the coding gain determines the 
number of users that can be supported by a base station. When 
this limit has been reached, a new user is denied access to the 
network and ripple effects of change is contained within the 
cell of the base station. 

In a wireless network without centralized control, there is no 
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Fig. 1. Illustration of the sand pile model. Each height plateau is a variable, 
and each variable can support up to z c sand grains, in this case 3. 



network entity to deny users access to the network resources. 
A simple example is Wi-Fi networks. Although too many Wi- 
Fi access points (APs) might be located in a given region 
for optimal operation, a person can set up a new Wi-Fi access 
point without any of the other users having a right to complain. 
For optimal operation after this new WIFI AP has been set 
up, the surrounding APs might have to change their transmit 
parameters such as power and frequency. 

B. An Example of a Self-organized Critical System 

A simple example of a dynamic system exhibiting SOC 
is the "sand pile" model JS). Assume we have N variables, 
where each variable z n is set to zero initially. Each variable 
is thought of as a discrete height plateau in a sand pile. A 
grain of sand is added to one variable at random. When z n is 
greater than a critical value z c , the grain tumbles to the next 
variable z n +i. This is illustrated in Fig. [TJ where z c — 3. In 
this example, if one grain is added to either z\ or z^ the grain 
tumbles down to z%, sliding either 1 or 2 plateaus down. In 
this system the length of a tumble, i.e. how many variables 
that surpass the critical value, also known as the avalanche 
size, behaves according to a power law. 

We can exchange the sand pile model with a set of wireless 
links, where we assume each link transmits with a power 
sufficient to support a given SINR. Then assume one link 
changes its power to support a new SINR requirement. If this 
change is large enough, "a grain tumbles" meaning some links 
update their power due to the change. These links are affected 
in the first "wave". If the accumulated change in power at 
these links exceeds a given threshold, a new set of links are 
affected and the avalanche continues to roll. 

II. System Model and Notation 

We consider a single frequency network where the decision 
variable at all users is power. We assume each user consists of 
a transmitter and a receiver, such that transmitter i transmits to 
receiver i. We denote the set of users as j\f, where N = \Af\. 
We assume each user i has a SINR requirement and 
transmits with a minimum power pi such that the SINR 
requirement is satisfied: 

l(d t i)pi 



where d^ is the distance between transmitter i and receiver 
j, l(dij) is the channel gain between i and j and No is the 
ambient noise variance. 

We use this allocation algorithm due to the fact that when 
a feasible allocation exist, this minimal power allocation will 
converge [lj. This way we are able to separate the domino 
effect from the cases when ripples of change occur due to the 
network not converging. 

Let pi be the transmit power of user i after the network 
has converged. Next, we assume a random user x updates its 
transmit power such that p' x = p x + A. Let p\ be the transmit 
power of user i after the network has converged again after 
user x's update. We are interested in characterizing the set 



A = {i\p'i ±Pi}. 



(2) 



Especially, we are interested in the cardinality of A. We can 
divide A into non-overlapping subsets A\, A2---, where the 
users belonging to Ai are the users directly affected by user 
x's update, i.e. are affected in the first round. The users 
belonging to A2 are the users not directly affected by user x's 
update, but due to the accumulated change in transmit power 
in the network from user x and the users in A\, changes its 
allocation, i.e. are affected in the second round. 

Characterizing A is difficult for multiple reasons. When 
simulating a network, multiple variables may be modeled 
as being drawn from different probability distributions. For 
instance, for a user i to belong to A\ we have that 
,1 



Pi 



Pi > Si 



(3) 



where Si is the threshold for user i to change its power. After 
user x updates its power, user i must change its power as 
follows to satisfy its SINR requirement 

1 _ Pi(N + lj + Al(d xi )) _ , R J{d x A 

Pi — u j \ — Pi ~T Pi^- 



l(du) 1,1 ' ^ l(d u ) ^ 

where Ii is the old interference at user i as given in the 
denominator in ([TJl. And thus user i changes its power if 

l(d x 
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All of these variables may be modeled according to some 
probability distribution. 

Throughout this paper we assume that the path loss function 
1(d) is a deterministic function of the distance d: 

l ( d ) = ^ ( 6 ) 
where a is the path loss exponent. 

A. Assumptions and Limitations 

We mentioned that the minimum power allocation scheme 
is used because when a feasible solution exists, this scheme 
will converge. Thus we are able to separate the domino 
effect from the cases when ripples of change occur due to 
the network not converging. This places constraints on the 
parameter values used in the simulation results throughout this 
paper. For instance, the density of a network seems to be an 




Fig. 2. Illustration of the array model. Distance between transmitter i and 
receiver i, da is the same for all i and dji = dj+i i+i = d«-_n_i. 

important parameter for the domino effect. It is reasonable to 
assume that dense networks will have larger domino effect than 
sparse networks, as users are closer to each other. However, 
we have not been able to simulate this aspect properly as 
with denser networks we need smaller SINR targets for the 
network to converge. For density issues, theoretical analysis 
seems actually more tractable. 

Two other aspects are related to the transmit power dis- 
tribution of the users. First, in practice each user has some 
maximum transmit power constraint. This will limit some 
users' ability to optimally adjust power after an update and 
also limit the domino effect. However, it is difficult to ana- 
lytically determine which users are operating near or at the 
maximum power and thus the analytical work in this paper 
is an upper bound on the domino effect. In the simulation 
results we enforce a maximum power constraint to obtain more 
realistic results. 

Secondly, according to (|4|)-(|5j, whether or not a user is 
affected by a change is determined by a difference in old and 
new transmit power. In the subsequent sections it will become 
apparent why this is done. However, from an engineering 
perspective it would be more useful to have such a criteria 
in dB. This again depends on the transmit power distribution 
of the users. E.g. two users may be transmitting at 0.01 W 
and 0.1 W. Both might have a difference of 0.05 between 
old transmit power and new transmit power after a change. 
However, for the first one this corresponds to a change of 7.8 
dB, for the other 0.2 dB change. 

III. Domino Effect in a Deterministic Model 

To illustrate the concept of the domino effect we start with 
a simple network model where transmitters and receivers are 
located along a line segment as illustrated in Fig. [2] We assume 
the number of transmitter and receiver pairs are infinite. At a 
time t transmitter updates its power by 1 unit. We assume 
that at time t < to the transmit power at each transmitter has 
converged such that changes in the network only occur due to 
the change at transmitter at time to. 

We let fti = 1 and da — d for all i and set A = 1. From 
Fig. [2] we have that 

dij+i = %dij (7) 

d ij+2 = yfiktfa - 2d 2 . (8) 

However, to simplify the expressions we assume dij = 
dij+i/Z = Gtj+2/3.... 

We now have that transmitter 1 will increase its transmit 
power if l(doi)/l(d) > 6. Transmitter 2 will increase its 
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Fig. 3. Affected number of users as function of cii and three values of <5 
(0.1 = -10 dB, 0.01 = -20 dB and 0.001 = -30 dB). The path loss 
exponent a = 3. 

transmit power if 

til a \ 1 ha iu a\ l ( d oi) . l(doi) 2 c 

(l(d 02 ) + l {dl2 ) — )/l{d) = — + — > S. 

We see that condition for whether transmitter n will increase 
its power becomes a sum depending on Z(doi) an d 1(d). To 
simplify notation, let a\ — l(doi)/l(d), which defines how 
closely the links are located. Assuming all transmitters up to 
n — 1 have increased their power, transmitter n will increase 
its power if 

n-l 

j=0 y ■>> 

For a given a\ and 6 we want to know if there exists a n such 
that a n < 5. This would mean that the change of transmitter 
affects a finite number of links, although this number may 
be large. 

Fig. [3] shows the number of affected users as a function of 
ai and 6 for a path loss exponent equal to 3. Decreasing 5 
leads more users to be affected by a change. We also see that 
for each S there exists a divergence value for a\ for which the 
number of affected users goes to infinity. 

IV. Poisson Distributed Network 

We now let the set of users be distributed over a 2- 
dimensional area according to a homogeneous Poisson Point 
Process (PPP) $. Specifically we let the transmitter locations 
be distributed according to a homogeneous PPP with density 
A. To avoid the frequency at which there is no feasible 
allocation, we let the receiver locations be dependent on 
the transmitter locations. Thus the receiver locations are also 
distributed according to a PPP, but these two distributions are 
not independent. 

A. Characterizing A\ 

To illustrate how changes in the network can affect different 
number of users we simplify the model by setting /?j = 1, 
da = 1/10 and Si = 0.01 for all i and set A = 1. 

The probability of user i belonging to A\ is thus the prob- 
ability that the path loss between transmitter x and receiver i 




Number oi Chan^i'i'l I 'sits 



Fig. 4. Analytical and simulated distributions of Ai referenced against 
simulated A with A = 4 X 10~ 4 , ft = 1, d„ = 1/10, a = 3 and Si = -20 
dB for all i and A = 1. 





Fig. 5. The plane divided into squares on the left hand side. On the right hand 
side, each square is assigned to an edge (bold lines) of a bond percolation 
model. 



is larger than 10 -5 . By assuming a — 3, it is the probability 
that transmitter x is closer than 46.4 meters to receiver i. A 
property of Slivnyak's theorem for PPPs is that the properties 
seen from a point x E R 2 are the same whether we condition 
on having a point of the PPP at x or not [9|. We therefore 
have 

j^l = $(46.4 2 7rA\{a;}) = $(46.4 2 ttA) (10) 
E[|^i|] = E[$(46.6 4 7rA\{a;})] = E[$(46.4 2 ttA)]. (11) 

Fig. |4] shows the analytical and simulated distribution of 
|^4i| for a density of nodes of 4 x 10 4 , referenced against the 
simulated distribution of \A\. As expected i?[|-4|] > 
and the probability of inflicting a large change in the network 
is larger for |.A| than for \A\\. 

B. Continuum Percolation and Infinite Propagation of Change 

From the estimate of \A\\ in the previous section we saw 
that E[|y^i |] is finite for all finite values of A. As expected, |.4| 
differed from |.4i| and we would like to characterize in 
some way. Unfortunately we are not able to provide any results 
on the distribution or moments of \A\, but we are able to 
provide results regarding the finiteness/infiniteness of \A\. To 
do this we use standard techniques from continuum percolation 
[9 1 to map the random Poisson network on a lattice and use 
results from bond percolation [10|. The result is as follows: 

Theorem 1. There exists a critical value for A, A c , for which 
values of A < A c the propagation of change is almost surely 
(a.s.) finite, and for values of X > A c the propagation of 
change is infinite a.s. 



Proof: Define r as 

r(£) = max{c? : /3£ 



Kd) 



>S} 



(12) 



for some variable £. We now say that two points i and j are 
connected if the distance between them is less or equal 
to r(£) for some £. We prove the theorem in two parts: first 
the existence of infinite propagation of change for large A and 
secondly the absence of infinite propagation of change for 
small A. 

1 ) Existence of infinite propagation of change for large X: 
Assume that the power update of all nodes except the initial 
node is S. If we can prove the existence of infinite propagation 
of change in this case, it must also hold for other update 
values as S is the lowest possible value. We now map the 
model onto a bond percolation model. The plane is divided into 
squares of size c = r(5) /2y/2 as shown in Fig. [5] Each square 
corresponds to a potential edge of the bond percolation lattice 
of Z 2 . The edge is added if at least one point of the Poisson 
process falls into this square. Each edge therefore exists with 
probability p = 1 — exp(— Ac 2 ), which is independent of all 
other edges. It is known that the bond percolation model of 
the infinite square lattice in Z 2 percolates if the edges between 
nearest neighboring points of the lattice exist with probability 
P > 1/2 |T0). Thus, if A > log2/c 2 , the edge percolation 
model contains an infinite cluster a.s. 

Now, if two edges are adjacent, it means that at least two 
points of the Poisson process are located in squares that share 
at least one corner, and therefore that the maximum distance 
between the two points is 2\p2c = 2\/2r(5) j2\[2 = r(S) and 
thus they can affect each other in the original model. As an 
infinite set of connected edges in the bond percolation model 
(i.e. infinite cluster) corresponds to an infinite set of connected 
points in the original model, there exists a path of infinite 
propagation of change a.s. if A > log2/c 2 . 

2) Absence of infinite propagation of change for small 
X: Let r be the distance r = r(A). Similarly, let r n 
be the distance r n = r(A„), where A„ will be defined 
below. Consider an initial node o that does an update A. 
The nodes initially affected by this change belong to the set 
A\. The expected number of nodes in A\ is Xnr^. Now, set 
A < 1/(7ITq) such that the expected number of affected nodes 
by the initial update is less than one. We now assume that 
only one node is affected in each "round", i.e. \A n \ — 
The maximum distance r n from a node in A n to a node in 
Ar, 



_i is r(A n ) and we have that A r 
A 



A„ = 2 Q ( 



n-l 



is given as 

5 



r-i 



+ S) (13) 



The factor 2 a comes from the fact that with a maximum radius 
r, the expected distance of a node within this area is r/2. 

'This does not affect the total number of affected users. As an example 
consider the case where there would be 2 affected users in one round, i.e. 
\A. n \ = 2. These two are affected independently of each other and we can 
divide this set into two sets .4* and A.„ . Now, all nodes that would be affected 
in An+l will still be affected as the accumulated change is the same. 



The question is now whether r m is finite or not. If r m is 
infinite then we will not be able to avoid percolation even 
for sufficiently small A, as long as A > 0. Since r n is just a 
product depending A n , r„ is finite if A„ is finite. 



The last term in ( 13 1 is a constant so we can drop this term. 
Also, for sufficiently large n the first term is zero. Hence, we 
focus on the sum Y^=i ^/(EjL; 1 r j) a as n — > oo. We also 
have that 
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when rj < rj+i,Vi > 1. Now, let 



i=2 



(14) 



(15) 



(16) 



The series £ converges if aj > i,Vi > 1 and a > 1, since 
Efcjv-'-A™ converges for any AT > when a > 1. dj > 
z, Vi > 1 holds whenever r2 > 2. 

Thus, as long as r2 > 2, r n < oo for all n and we can 
choose A small enough so that the average number of affected 
nodes in each round is always less than 1. This process is 
equivalent to a Galton- Watson process and it is known that 
when the average number of children per individual in a 
Galton-Watson process is smaller than one, the process dies 
out with probability 1 ifTTIl . 

Remark: The inequality between ( fT~4-| > and ( 15 1 holds as long 
as j-j < ri + i,Vi > 1. If this is not the case for some i, then 
we have a maximum value of r, for some i < oo. Since i 
must be less than infinity, this maximum must also be less 
than infinity and we can choose A small enough so that the 
average number of affected nodes from this distance is less 
than one. Thus the result holds. ■ 

Although we have proved that a critical density exists, we 
have not provided any value for it. From the proofs of both 
existence and absence of percolation we also see that such a 
critical density depends on the other system parameters fa, da 
and S. For fixed fa and da, the critical density increases with 
S. 

C. Simulation Results and Self-organized Criticality 

As we have not been able to characterize \A\ analytically, 
we investigate the distribution of \A\ through simulations 
in this section. However, we do have some expectations of 
the distribution of \A\ depending on the distributions of the 
individual |^|s. If |A| << for i > 1, we expect \A\ to 
be almost Poisson distributed since |.Ai| is Poisson distributed. 
If E[|A|] > it is less clear how \ A\ should be distributed. It 
is well known that the sum of independent Poisson variables 
are still Poisson. But in this case j^l is not independent of 
| Ai \, since if a realization of j^4_x | is small, we also expect \A2 1 
to be small. Whether or not the sum is Poisson depends on the 
conditional probability function [12|. Through the simulation 
results, it seems like the resultant distribution is still Poisson. 
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Fig. 6. Distribution of |_4| for different values of <5 and A = 4 X 10 4 . 




Fig. 7. Distribution of the different \Ai\sforS = -30 dB and A = 4xl0~ 



We have used the same parameters as in Section IV-A 
namely A = 4 x 1CT 4 , fa = 1, d u = 1 0, A = 1, a = 3, 
No = 10~ 8 . As mentioned in Section II-A we did not assume 



a maximum power constraint in the analytical work in the 
previous sections. In the simulations we assume a maximum 
transmit power of 1 W, and thus even though a transmitter 
should update its power according to |5]), if the maximum 
power will be surpassed by such a change the transmitter will 
not update its power. 

In Fig.|6]the distribution of \A\ is plotted for different values 
of S. For high values of S, \Ai\ << \Ai\, for i > 1 and 
the plots for i5 = —10 dB and —20 dB are almost Poisson 
distributed. 

In Fig. [Tjthe distribution of the different |*Ai|s for S = —30 
dB and A = 4x 10~ 4 is given. We see that all the |A;|s seem 
to follow a Poisson distribution with different expected values. 

In Fig. [8] we have plotted which Ai a fraction of A belongs 
to. We see that this distribution follows a power law for small 
5. In relation to self-organized criticality, we see that the size 
of avalanches, where avalanche size is given as the maximum 
i for which \Ai\ is nonzero, follows the same characteristics 
of a system in a critical state. However, for 6 = — 10 dB, the 
curve does not follow a power law distribution and we expect 
this state to be non-critical. 

From the results presented above, we see that a large number 
of users can be affected by a change in a manner similar to a 
domino effect. The number of users that are affected depends 
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Fig. 8. Fraction of |.4| in the different |*4i|s in a loglog plot for different 
values of 8 and A = 4 X 10 -4 . 




Fig. 9. Average power increase in dB for different values of <5 with A 
4 X 1CT 4 . 



Through simulation results we showed how the number of 
affected nodes follows different Poisson distributions in a 
network distributed according to a Poisson Point process, that 
depends on the threshold for when a node is considered to be 
affected. While the domino effect is a decreasing function of 
threshold, one cannot simply set a high threshold as this leads 
to significant increase in power and degradation of system 
performance. Thus the optimal threshold is a trade-off be- 
tween having high system performance while still maintaining 
resilience to the domino effect. 

While this paper has highlighted some of the issues regard- 
ing change in decentralized wireless network, there are still 
important open issues. One is to give upper and lower bounds 
on the critical density for percolation. Another is to give a 
characterization of the trade-off one gets from limiting the 
domino effect while maximizing system performance. 

Lastly, this paper only considered power control in a Poisson 
distributed network. We did this because a simple power 
control algorithm is guaranteed to converge when a feasible 
solution exists and hence we were able to separate change 
due to convergence and change due to the domino effect. An 
equally important aspect is how frequency allocation affects 
the domino effect. As frequency is a discrete variable it might 
in some cases decrease the propagation of change as it can be 
seen as a thinning of the Poisson process. On the other hand, 
a slight change in power at one transmitter might result in a 
frequency change at another which can escalate the domino 
effect. 



on the threshold <5. For analytical purposes it was convenient to 
define the condition for when a user has to update its power |5]) 
on a difference between old power and new necessary power. 
However, from an engineering aspect it would be better to have 
this in terms of dB. In Fig. [9] the average power increase in 
dB is shown for different values of S for the same simulation 
parameters used in results Fig. 6|8 As we see, even a modest 
threshold of S = —30 dB, leads to an average power increase 
over all users of 18 dB. This point out the difficulty in choosing 
a suitable threshold, as a high threshold leads to a low domino 
effect, but leads to an increased power which degrades the 
performance. 

V. Conclusion and Future Work 

In this paper we have investigated the issue of propagation 
of change in decentralized wireless networks, an issue largely 
lacking in the literature. This manifests itself as a domino 
effect in the network can arise due to new users entering 
the network or old ones leaving or updating their resource 
allocation. As most proposed resource allocation algorithms 
are dynamic and adaptive, they change according to the current 
state of the network in order to utilize the resources in an 
optimal manner. However this leads to ripples of change 
propagating through the network. In this paper we showed 
that there exists a critical network density, for which a larger 
density will propagate change infinitely almost surely, while 
smaller densities will propagate change finitely almost surely. 
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